Tunneling probability for the birth of an asymptotically DeSitter universe. 
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In the present work, we quantize a closed Friedmann-Robertson- Walker model in the 
^Hy presence of a positive cosmological constant and radiation. It gives rise to a Wheeler-DeWitt 

equation for the scale factor which has the form of a Schrodinger equation for a potential 



with a barrier. We solve it numerically and determine the tunneling probability for the 
birth of a asymptotically DeSitter, inflationary universe, initially, as a function of the mean 
energy of the initial wave-function. Then, we verify that the tunneling probability increases 
with the cosmological constant, for a fixed value of the mean energy of the initial wave- 
function. Our treatment of the problem is more general than previous ones, based on the 
WKB approximation. That is the case because we take into account the fact that the scale 
factor (a) cannot be smaller than zero. It means that, one has to introduce an infinity 
potential wall at a = 0, which forces any wave-packet to be zero there. That condition 
introduces new results, in comparison with previous works. 
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I. INTRODUCTION 

Since the pioneering work in quantum cosmology due to DeWitt [1], many physicists have worked 
in this theory. The main motivation behind quantum cosmology is a consistent explanation for 
the origin of our Universe. So far, the most appealing explanation is the spontaneous creation 
from nothing In that picture for the origin of the Universe, the Universe is a 

quantum mechanical system with zero size. There is a potential barrier that the Universe may 
tunnel with a well-defined, non-zero probability. If the Universe actually tunnels, it emerges to 
the right of the barrier with a definite size. The application of the creation from nothing idea in 
minisuperspace models has led to several important results. The wave-function of the Universe 
satisfies the Wheeler-DeWitt equation [l[ |J. Therefore, one needs to specify boundary conditions 
in order to solve that equation and find a unique and well-defined wave function. The motivation 
to obtain a wave-function that represents the creation from nothing has led to the introduction of 
at least three proposals for the boundary conditions for the wave- function of the Universe [3]. The 
inflationary period of the Universe appears very naturally from the creation from nothing idea. 
That is the case because most of the minisuperspace models considered so far have a potential 
that decreases, without a limit, to the right of the barrier. It gives rise to a period of unbounded 
expansion which is interpreted as the inflationary period of the Universe Q]. Also, it was shown by 
several authors that an open inflationary universe may be created from nothing, in theories of a 



single scalar field for generic potentials 
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10, 



111 ]. Another important issue is the particle content 
in the Universe originated during the creation from nothing process H . . 

In the present work, we would like to explicitly compute the quantum mechanical probability 
that the universe tunnels through a potential barrier and initiates its classical evolution. That 
probability is the tunneling probability (TP) and the particular model we consider, here, is a closed 
Friedmann-Robertson- Walker (FRW) model in the presence of a positive cosmological constant and 
radiation. The radiation is treated by means of the variational formalism developed by Schutz []jj]. 
That model has already been treated quantum mechanically using the ADM formalism and the 
Dirac quantization for constrained systems [?], Q, Q]- The wave- function, for that model, was 



calculated in the WKB approximation. Here, we compute the wave-function and TP exactly, 
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without any approximation. It will be done by means of a numerical calculation. In particular, 
our treatment of the problem is more general than previous ones, because we take into account 
the fact that the scale factor (a) cannot be smaller than zero. It means that, one has to introduce 
an infinity potential wall at a = 0, which forces any wave-packet to be zero there. As we shall see, 
that condition introduces new results, in comparison with previous works. This model has two 
free parameters: the radiation energy and the cosmological constant. Therefore, we will obtain 
the TP as a function of those two parameters. One of the main motivations of any quantum 
cosmological model is to fix initial conditions for the classical evolution of our Universe [2]. Here, 
for the present model, we would gain information on what is the most probable amount of radiation 
in the initial evolution of the classical universe and the most probable value of the cosmological 
constant. Another motivation of the present work, is trying to contribute to a long standing debate 
about which is the most appropriate set of initial conditions for the wave-function of the Universe. 
The most well-known proposals for the set of initial conditions are the tunneling one, due to A. 
Vilenkin [3], and the no-boundary one, due to J. B. Hartle and S. W. Hawking [Jj. The application 
of those proposals for simple models showed that they give some different predictions for the initial 



evolution of the Universe H,MLI13, 15, 
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181 ] . One of such predictions, which we shall explore 



here, is the initial energy of the Universe right after its nucleation. The tunneling wave function 
predicts that the Universe must nucleate with the largest possible vacuum energy whereas the 
no-boundary wave function predicts just the opposite [7]. In terms of our results, if one assumes 
that the cosmological constant describes a vacuum energy, it is important to see if TP increases or 
decreases with the cosmological constant, for fixed radiation energy. The first behavior favors the 
tunneling wave function and the latter favors the no-boundary wave function. 

In the next Section, we describe the classical dynamics of the present cosmological model. We 
write the super-hamiltonian constraint and the Hamilton's equations. We solve the Hamilton's 
equations and find the general solution of the model. Then, we qualitatively describe all possible 
classical evolutions. In Section HH] we canonically quantize the model and obtain the corresponding 
Wheeler-DeWitt equation. We solve it, numerically, for particular values of the radiation energy 
and the cosmological constant. We show the square modulus of the wave-function of the universe 
as a function of the scalar factor. The tunneling process can be readily seen from that figure. 
The Section IIVI is divided in three subsections with the main results of the paper. In the first 
subsection IIV Al we start introducing the tunneling probability, then we evaluate its dependence 
on the radiation energy. We obtain that, the TP increases with the radiation energy for a fixed 
cosmological constant. Therefore, it is more probable that the classical evolution should start 
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with the greatest possible value for the radiation energy. In the following subsection IIVBI we 
give a detailed comparison between the exact TP, computed in the previous subsection, and the 
corresponding WKB tunneling probability. Here, we show how the presence of an infinity potential 
wall at a = may lead to a difference between our results and previous ones, based on the WKB 
approximation. In the final subsection IIV CI of this section, we evaluate the dependence of the 
exact TP with the cosmological constant. We obtain that, the TP increases with the cosmological 
constant for a fixed radiation energy. Therefore, it is more probable that the classical evolution 
should start with the greatest possible value for the cosmological constant. This behavior of TP 
also favors the tunneling wave function. Finally, in Section |V] we summarize the main points and 
results of our paper. 



II. THE CLASSICAL MODEL 



The Friedmann-Robertson- Walker cosmological models are characterized by the scale factor 
a{t) and have the following line element, 

ds 2 = -N 2 (t)dt 2 + a\t) + , (1) 

where dQ, 2 is the line element of the two-dimensional sphere with unitary radius, N(t) is the 
lapse function and k gives the type of constant curvature of the spatial sections. Here, we are 
considering the case with positive curvature k = 1 and we are using the natural unit system, where 
h = c = G = 1. The matter content of the model is represented by a perfect fluid with four- velocity 
U 11 = 5q in the comoving coordinate system used, plus a positive cosmological constant. The total 
energy-momentum tensor is given by, 



T (l>v = {p + p)U ll U v -pg^ v -Kg lhV , (2) 

where p and p are the energy density and pressure of the fluid, respectively. Here, we assume that 
p = p/3, which is the equation of state for radiation. This choice may be considered as a first 
approximation to treat the matter content of the early Universe and it was made as a matter of 
simplicity. It is clear that a more complete treatment should describe the radiation, present in the 
primordial Universe, in terms of the electromagnetic field. 

Einstein's equations for the metric (pQ) and the energy momentum tensor ([2]) are equivalent to 
the Hamilton equations generated by the super-hamiltonian constraint 
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W = __^L_ 3 a + Aa 3 + ^, 
12a a 



(3) 



where p a and pt are the momenta canonically conjugated to a and T the latter being the canonical 
variable associated to the fluid [lj3]. The total Hamiltonian is given by NTL and we shall work in 
the conformal gauge, where N = a. 

The classical dynamics is governed by the Hamilton equations, derived from the total Hamilto- 
nian N7i, namely 



d(NH) 



Pa 



Pa 



-^)=6a-4Aa3 
oa 



Pt 



d(NH) 
c>Pt 

d(NH) 
dT 



1, 



(4) 



0. 



Where the dot means derivative with respect to the conformal time r = Nt. We also have the 
constraint equation TL = 0. We have the following solutions for the system @: 



T(r) = T + ci,a(r) 



-.sn 



18 + 6 ^9- 12A/5 (r - r ) 
6 



(5) 



'3+ y/9- V2A/3 

where ci, /? and To are integration constants, sn is the Jacobi's elliptic sine |2U] of modulus rr givcu 
by 



v 7 ^ /-2/3A + 3- V9 - 12 A ^ 



(6) 



2 V A/3 

In the case studied here A > 0, the radiation energy (5 can take values in the domain, (3 < 3/(4A). 
If one substitutes values of (3 such that (3 < 3/(4A) in Eqs. © and ([6]), the scale factor, starting 
from zero, expands to a maximum size and then recollapse. On the other hand, if (3 = 3/(4A), the 
scalar factor initially decelerates and then, enter the regime of unbounded expansion. 
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III. THE QUANTUM MODEL 

We wish to quantize the model following the Dirac formalism for quantizing constrained systems 



2l\ . First we introduce a wave- function which is a function of the canonical variables a and T, 



V = *(a,T). (7) 

Then, we impose the appropriate commutators between the operators a and T and their conjugate 
momenta P a and Pt- Working in the Schrodinger picture, the operators a and T are simply mul- 
tiplication operators, while their conjugate momenta are represented by the differential operators 

d d 

Pa - , PT - -i^ ■ (8) 

Finally, we demand that the operator corresponding to N7i annihilate the wave- function *S>, 
which leads to Wheeler-DeWitt equation 

3 ° 2 + Afl4 ) * (a ' r) = _< J^ (a ' r) ' (9) 

where the new variable r = — T has been introduced. 

The operator NTL is self-adjoint [22j with respect to the internal product, 



POO 

(¥,$)=/ dam{a,TY^(a,T), (10) 
Jo 

if the wave functions are restricted to the set of those satisfying either ^(0, r) = or ^'(0, r) = 0, 
where the prime / means the partial derivative with respect to a. Here, we consider wave functions 
satisfying the former type of boundary condition and we also demand that they vanish when a 
goes to oo. 

The Wheeler-DeWitt equation Q is a Schrodinger equation for a potential with a barrier. 
We solve it numerically using a finite difference procedure based on the Crank-Nicholson method 
231 ]. 24 ] and implemented in the program GNU-OCTAVE. Our choice of the Crank-Nicholson 
method is based on its recognized stability. The norm conservation is commonly used as a criterion 
to evaluate the reliability of the numerical calculations of the time evolution of wave functions. 



In References 25( and 2g], this criterion is used to show analytically that the Crank-Nicholson 
method is unconditionally stable. Here, in order to evaluate the reliability of our algorithm, we 
have numerically calculated the norm of the wave packet for different times. The results thus 
obtained show that the norm is preserved. 
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In fact, numerically one can only treat the tunneling from, something process, where one gives a 
initial wave function with a given mean energy, very concentrated in a region next to a = 0. That 
initial condition fixes an energy for the radiation and the initial region where a may take values. 
Our choice for the initial wave function is the following normalized gaussian, 



where E is the radiation energy. ^f(a, 0) is normalized by demanding that the integral of \^(a, 0)| 2 
from to oo be equal to one and its mean energy be E. After one gives the initial wave function, 
one leaves it evolve following the appropriate Schrodinger equation until it reaches infinity in the 
a direction. Numerically, one has to fix the infinity at a finite value. In the present case we fix 
Umax = 30 as the infinity in the a direction. The general behavior of the solutions, when E is 
smaller than the maximum value of the potential barrier, is an everywhere well-defined, finite, 
normalized wave packet. Even in the limit when the scale factor goes to zero. For small values of 
a the wave packet have great amplitudes and oscillates rapidly due to the interaction between the 
incident and reflected components. The transmitted component is an oscillatory wave packet that 
moves to the right and has a decreasing amplitude which goes to zero in the limit when a goes 
to infinity. As an example, we solve eq. ([9]) with A = 0.0121. For this choice of A the potential 
barrier has its maximum value equal to 185.95. In order to see the tunneling process, we choose 
E = 185 for the initial wave function eq. (jlip. For that energy, we compute the points where it 
meets the potential barrier, the left (c% p ) and right (a r t p ) turning points. They are, au v = 10.7287 
and a r t p = 11.5252. In figure [TJ we show |^(a, t max )\ 2 for the values of A and E, given above, at 
the moment (t max ) when ^ reaches infinity. For more data on this particular case see table Hill in 
the appendix. It is important to mention that the choice of the numerical values for A and E above 
and in the following examples, in the next section, were made simply for a better visualization of 
the different properties of the system. 
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IV. TUNNELING PROBABILITIES 



A. Tunneling probability as a function of E 



We compute the tunneling probability as the probability to find the scale factor of the universe 
to the right of the potential barrier. In the present situation, this definition is given by the following 
expression, 



8 



0.4 



P 



0.2 



0.3 



0.1 




5 



10 



15 



20 



25 



30 



a 



FIG. 1: \^f(a, t max ) 2 = p, for A = 0.0121, E = 185 at the moment t max when $ reaches infinity, located at a = 30. 



where, as we have mentioned above, numerically oo has to be fixed to a maximum value of a. Here, 
we are working with dmax ~ 

30. 

Since, by normalization, the denominator of Eq. (|12p is equal to the identity, TPi n t is effectively 
given by the numerator of Eq. (|12|) . We consider initially the dependence of TP on the energy 
E. Therefore, we compute TPi n t for many different values of E for a fixed A. For all cases, we 
consider the situation where E is smaller than the maximum value of the potential barrier. From 
that numerical study we conclude that the tunneling probability grows with E for a fixed A. As 
an example, we consider 47 values of the radiation energy for a fixed A = 0.01. For this choice 
of A the potential barrier has its maximum value equal to 225. In order to study the tunneling 
process, we fix the mean energies of the various ^(a, 0)'s eq. (|lip to be smaller than that value. 
In table HU in the appendix, we can see, among other quantities, the different values of the energy 
E, TPi n t, OLitp an d cirtp for each energy. In figure [21 we see the tunneling probability as functions 
of E, for this particular example. Due to the small values of some TP's, we plot the logarithms of 
the TP's against E. 

Since TP grows with E it is more likely for the universe, described by the present model, to 
nucleate with the highest possible radiation energy. Therefore, it is more probable that the classical 
evolution should start with the greatest possible value for the radiation energy. 



TPint 



IZ P \^(a,tmax)\ 2 da 




(12) 



Jq°° \^(a,tmax)\ 2 da 
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FIG. 2: \ogTPint for different radiation energies (E) for a fixed A = 0.01. 



B. Exact tunneling probability versus WKB tunneling probability 



Let us, now, compare the exact tunneling probability represented by TPi n t eq. (|12p with 
the approximated WKB tunneling probability (TPwkb)- The TPy/KB is defined by the ratio 
between the square modulus of the transmitted amplitude of the WKB wave-function and the 
square modulus of the incident amplitude of the WKB wave- function [3], [28]. For the present 



situation, it is given by the following expression 



28] 



TPwKB = ~2 > (13) 



(20 



+ 



where, 

= exp ( da^Jl2(3a 2 - Aa 4 - E^J . (14) 

It is important to note that the TPwkb eq. (fT3|). was computed considering that the incident 
wave (\I//) reaches the potential barrier at au p . Then, part of *$>i is transmitted to 00 (^>t) and 
part is reflected to — 00 (^r). In the present problem, we have an infinity potential wall at a = 
because the scale factor cannot be smaller than zero. It means that cannot go to —00, as 
was assumed in order to compute the TPwkb- Instead, will reach the infinity potential wall 
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at a = and will be entirely reflected back toward the potential barrier, giving rise to a new 
incident wave ({^r)i). The new incident wave r)i reaches the potential barrier at au p and is 
divided in two components. A reflected component which moves toward the infinity potential wall 
at a = (((^r)i)r) and a transmitted component which moves toward oo ({(^?r)i)t)- ((^r)i)t 
will contribute a new amount to the already existing TPi n t due to (^t). On the other hand, 
((^r)i)r, after being reflected at a = 0, gives rise to a new incident wave which will contribute a 
further amount to the already existing TPi n f If we let our system evolve for a long period of time, 
TPi n t will get many such contributions from the different reflected components. The only way it 
makes sense comparing TP{ n t with TPwkb is when we let the system evolve for a period of time 
(At) during which cannot be reflected at a = and come back to reach the potential barrier. 
It is clear by the shape of our potential that the greater the mean energy E of the wave-packet 
(jlip . the greater is (At). As an example, in Table HI we show a comparison between TPi n t and 
TPwkb for different values of E and At for the case with A = 0.01. We can see, clearly, that both 
tunneling probabilities coincide if we consider the appropriate At, for each E. 



E 


TPint 


TPwkb 


At 


80 


6.0648 x 10~ 302 


1.1845 x 10~ 303 


13 


100 


1.5887 x 10" 259 


3.7271 x 10" 258 


18.5 


130 


6.4375 x lO" 194 


9.8051 x 10" 193 


30 


160 


6.9194 x 10" 130 


5.7774 x 10~ 130 


45.5 


175 


6.5061 x 10" 100 


2.1361 x 10"" 


54.5 


190 


5.3119 x 10~ 69 


2.6372 x 10" 69 


65.5 


200 


2.2682 x 10" 49 


1.7295 x 10" 49 


73.5 


215 


5.4531 x 10~ 20 


4.1983 x 10" 20 


88 


219 


5.2168 x 10" 12 


2.4754 x 10" 12 


93 


223 


7.0045 x 10 -04 


1.3731 x 10" 04 


100 



TABLE I: A comparison between TPi nt and TPwkb for 10 different values of E with its associated integration 
time At for the case with A = 0.01. 

In order to have an idea on how the TPi n t may differ from the TPwkb, we let the initial wave- 
packet (jlip . with different mean energies, evolve during the same time interval At. We consider 
the example given in the previous subsection IIV Al with a common time interval of 100. Observing 
Table HI we see that this amount of At guarantees that ^r of the wave-packet with mean energy 
223 does not contribute to the TP{ n t. Therefore, we may expect that to all wave-packets with mean 
energies smaller than 223, TPi n t will be greater than TPwkb- We show this comparison in table 
ILTl in the appendix, where we have an entry for TPwkb- It means that, we computed the values 
of the TPwkbs for each energy used to compute the TPj^s, in the case where A = 0.01. In figure 
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[3l we show, graphically, that comparison between both tunneling probabilities as functions of E, 
with At = 100 for all values of E. Due to the small values of some TP's, we plot the logarithms 
of the TP's against E. 
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FIG. 3: Comparison between logTPwKB (dots) and \ogTPi nt (line) for different radiation energies (E) for a fixed 
A = 0.01 . The integration time At is equal to 100 for all values of E. 

As we can see from figure El for this choice of At the tunneling probabilities disagree for most 
values of E. They only agree for values of E very close to the top of the potential barrier. There, 
because the values of E are similar to 223, At is almost sufficient to guarantee that *$>r of each 
wave-packet does not contribute to the TPi n t- 

As we have mentioned above, numerically one can only treat the tunneling from something 
process, where one gives a initial wave function with a given mean energy, very concentrated in a 
region next to a = 0. Then, if we take E = the TPi n t will be zero. On the other hand, we may 
have an idea how TPi n t behaves near E = from figure [21 From table HU in the appendix, one can 
see that TP WKB = 7.0246 x 10~ 522 when E = 0. 

Finally, we may compute the time (r) the universe would take, for each energy, to nucleate. 
In order to understand the meaning of r, consider a photon that composes the radiation which is 
initially confined to the left of the potential barrier. Then, compute the emission probability of 
that photon, per unit of time. We may invert it to obtain r, the time the photon would take to 
escape the potential barrier. If we consider r as the time it takes for the most part of the photons 
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to escape the barrier, we obtain the time the universe would take, for each energy, to appear at 



the right of the barrier. In the present situation, r is given by the following expression 



t = 2ai tp —— (15) 

From table HU in the appendix, we may see the values of r, for each energy. It is clear by the 
results that, the smaller the energy E the longer it will take for the universe to nucleate. 

C. Tunneling probability as a function of A 

We would like to study, now, how the tunneling probability depends on the cosmological con- 
stant. In order to do that, we must fix an initial energy E for the radiation and compute the 
TPi n t for various values of A. For all cases, we consider the situation where E is smaller than the 
maximum value of the potential barrier. From that numerical study we conclude that the tunneling 
probability grows with A for a fixed E. As an example, we choose E = 185 and 22 different values 
of A, such that, the maximum energy of the potential barrier (PE max ), for each A, is greater than 
185. The values of A, TPi nt , PE max , r, au p and a r t p are given in table in the appendix. With 
those values, we construct the curve TPi n t versus A, shown in figure 01 Due to the small values of 
some TP's, we plot the logarithms of the TP's against A. 

Since TP grows with A it is more likely for the universe, described by the present model, to 
nucleate with the highest possible cosmological constant. Therefore, it is more probable that the 
classical evolution should start with the greatest possible value for the cosmological constant. Also, 
if we assume that A describes a vacuum energy, this result is qualitatively in accordance with the 
prediction of the tunneling wave function due to A. Vilenkin [3]. 

V. CONCLUSIONS. 

In the present work, we canonically quantized a closed Friedmann-Robertson- Walker (FRW) 
model in the presence of a positive cosmological constant and radiation. The radiation was treated 
by means of the variational formalism developed by Schutz [ijj]. The appropriate Wheeler-DeWitt 
equation for the scale factor has the form of a Schrodinger equation for a potential with a barrier. 
We solved it, numerically, and determined the tunneling probability for the birth of an asymptot- 
ically DeSitter, inflationary universe, as a function of the radiation energy and the cosmological 
constant. In particular, our treatment of the problem is more general than previous ones, because 
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FIG. 4: log TP in t for 22 different values of A for a fixed E = 185. 



we took into account the fact that the scale factor (a) cannot be smaller than zero. It means 
that, one has to introduce an infinity potential wall at a = 0, which forces any wave-packet to 
be zero there. As we saw, that condition introduced new results, in comparison with previous 
works. Then, we verified that the tunneling probability increases with the radiation energy, for 
a fixed cosmological constant. Therefore, it is more probable that the classical evolution should 
start with the greatest possible value for the radiation energy. We also gave a detailed comparison 
between the exact TP, computed here, and the corresponding WKB tunneling probability. Finally, 
we evaluated the dependence of the exact TP with the cosmological constant. We obtained that, 
the TP increases with the cosmological constant for a fixed radiation energy. Therefore, it is more 
probable that the classical evolution should start with the greatest possible value for the cosmo- 
logical constant. Also, if one assumes that the cosmological constant describes a vacuum energy, 
the latter result seems to be in accordance with the predictions of the tunneling wave function 3J]. 
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Encrgia 


TPint 


TPwKB 


T 


altp 


artp 


0.0000 


0.0000 


7.0246 x 10~ 522 


oo 


0.0000 


17.3205 


1.0000 


2.5795 x 10~ 67 


2.7574 x 10~ 517 


4.4790 x 10+ 66 


0.5777 


17.3109 


2.0000 


3.9975 x 10- 64 


2.7181 x lO" 513 


4.0896 x 10+ 63 


0.8174 


17.3012 


3.0000 


4.8040 x 10~ 61 


1.5939 x lO" 509 


4.1702 x 10+ 60 


1.0017 


17.2915 


4.0000 


2.6388 x 10" 59 


6.6774 x 10" 506 


8.7714 x 10+ 58 


1.1573 


17.2818 


5.0000 


5.7738 x 10~ 57 


2.1799 x 10~ 502 


4.4844 x 10+ 56 


1.2946 


17.2721 


6.0000 


9.3459 x 10~ 56 


5.8369 x 10~ 499 


3.0366 x 10+ 55 


1.4190 


17.2623 


7.0000 


6.9178 x 10~ 55 


1.3258 x 10~ 495 


4.4337 x 10+ 54 


1.5335 


17.2525 


8.0000 


7.0061 x 10" 56 


2.6169 x 10" 492 


4.6827 x 10+ 55 


1.6404 


17.2427 


9.0000 


6.3878 x 10" 53 


4.5691 x l(r 489 


5.4506 x 10+ 52 


1.7409 


17.2328 


10.0000 


3.9939 x 10~ 51 


7.1563 x 10" 486 


9.1944 x 10+ 50 


1.8361 


17.2229 


15.0000 


1.3310 x 10~ 46 


1.8319 x 10~ 470 


3.3888 x 10+ 46 


2.2553 


17.1731 


20.0000 


3.3918 x 10~ 44 


9.0816 x 10- 456 


1.5401 x 10+ 44 


2.6119 


17.1224 


30.0000 


1.4814 x 10~ 41 


7.5933 x 10~ 428 


4.3450 x 10+ 41 


3.2183 


17.0189 


40.0000 


5.8991 x 10~ 40 


2.5466 x lO" 401 


1.2679 x 10+ 40 


3.7397 


16.9120 


50.0000 


9.8017 x 10~ 39 


8.0358 x 10" 376 


8.5875 x 10+ 38 


4.2086 


16.8014 


60.0000 


1.1252 x 10~ 37 


3.9314 x 10~ 351 


8.2507 x 10+ 37 


4.6419 


16.6869 


70.0000 


1.1121 x 10" 36 


4.1409 x 10~ 327 


9.0821 x 10+ 36 


5.0499 


16.5680 


80.0000 


1.0627 x 10" 35 


1.1845 x 10~ 303 


1.0236 x 10+ 36 


5.4391 


16.4443 


90.0000 


1.0557 x 10" 34 


1.0939 x 10~ 280 


1.1016 x 10+ 35 


5.8147 


16.3153 


100.0000 


1.1488 x 10~ 33 


3.7271 x 10~ 258 


1.0760 x 10+ 34 


6.1803 


16.1803 


110.0000 


1.4319 x 10" 32 


5.2113 x 10~ 236 


9.1338 x 10+ 32 


6.5393 


16.0386 


120.0000 


2.1333 x 10~ 31 


3.2602 x 10~ 213 


6.4634 x 10+ 31 


6.8942 


15.8893 


130.0000 


3.9754 x 10~ 30 


9.8051 x 10~ 193 


3.6464 x 10+ 30 


7.2479 


15.7311 


140.0000 


9.7584 x 10" 29 


1.5060 x 10~ 171 


1.5582 x 10+ 29 


7.6029 


15.5626 


150.0000 


3.3597 x 10" 27 


1.2439 x lO" 150 


4.7398 x 10+ 27 


7.9623 


15.3819 


160.0000 


1.7562 x 10~ 25 


5.7774 x 10" 130 


9.4858 x 10+ 25 


8.3293 


15.1863 


165.0000 


1.5321 x 10~ 24 


1.0157 x 10~ 119 


1.1118 x 10+ 25 


8.5171 


15.0818 


170.0000 


1.5472 x 10~ 23 


1.5685 x 10" 109 


1.1259 x 10+ 24 


8.7085 


14.9720 


175.0000 


1.8431 x 10" 22 


2.1361 x 10~ 99 


9.6624 x 10+ 22 


8.9045 


14.8563 


180.0000 


2.6520 x 10" 21 


2.5758 x 10~ 89 


6.8673 x 10+ 21 


9.1059 


14.7337 


185.0000 


4.7418 x 10- 20 


2.7600 x 10" 79 


3.9286 x 10+ 20 


9.3142 


14.6029 


190.0000 


1.0919 x 10~ 18 


2.6372 x 10~ 69 


1.7457 x 10+ 19 


9.5310 


14.4624 


195.0000 


3.3916 x 10~ 17 


2.2544 X 10~ 59 


5.7545 x 10+ 17 


9.7585 


14.3099 


200.0000 


1.5114 x 10~ 15 


1.7295 x 10" 49 


1.3233 x 10+ 16 


10.0000 


14.1421 


205.0000 


1.0542 x 10" 13 


1.1943 x 10" 39 


1.9466 x 10+ 14 


10.2605 


13.9543 


210.0000 


1.3129 x lO" 11 


7.4432 x 10~ 30 


1.6069 x 10+ 12 


10.5485 


13.7379 


215.0000 


3.6494 x 10~ 09 


4.1983 x lO" 20 


5.9628 x 10+ 09 


10.8801 


13.4767 


216.0000 


1.2796 x 10~ 08 


3.7003 x 10" 18 


1.7121 x 10+ 09 


10.9545 


13.4164 


217.0000 


4.7368 x lO" 08 


3.2487 x 10~ 16 


4.6582 x 10+ 08 


11.0325 


13.3523 


218.0000 


1.8642 x lO"" 7 


2.8413 x 10~ 14 


1.1926 x 10+ 08 


11.1150 


13.2837 


219.0000 


7.8683 x 10" 07 


2.4754 x 10~ 12 


2.8476 x 10+ 07 


11.2029 


13.2097 


220.0000 


3.6052 x 10-° 6 


2.1485 x lO" 10 


6.2674 x 10+ 06 


11.2978 


13.1286 


221.0000 


1.8228 x 10~ 05 


1.8577 x 10~ 8 


1.2512 x 10+ 06 


11.4018 


13.0384 


222.0000 


1.0419 x lO" 04 


1.6002 x 10~ 6 


2.2110 x 10+ 05 


11.5187 


12.9352 


223.0000 


7.0045 x 10" 04 


1.3731 x 10~ 4 


3.3281 x 10+ 04 


11.6558 


12.8118 


224.0000 


5.9816 x 10" 03 


1.1671 x 10" 2 


3.9562 x 10+ 03 


11.8322 


12.6491 



TABLE II: The computed values of TP in t, TPwkb, t, au p 



and a r t P for 47 different values of E when A = 0.01. 
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A 


TPint 


PE 


r 


a it P 


CLrtp 


0.0100 


A.7AA9 x 10~ 20 


2.2500 x 10+ 02 


3.9260 x 10+ 20 


9.3142 


14.6029 


0.0101 


1.9386 x 10~ 17 


2.2277 x 10+ 02 


9.6427 x 10+ 17 


9.3467 


14.4800 


0.0102 


5.5162 x 10~ 15 


2.2059 x 10+ 02 


3.4010 x 10+ 15 


9.3803 


14.3571 


0.0103 


1.0768 x 10" 12 


2.1845 x 10+ 02 


1.7487 x 10+ 13 


9.4152 


14.2343 


0.0104 


1.4239 x 10- 10 


2.1635 x 10+ 02 


1.3276 x 10+ 11 


9.4515 


14.1114 


0.0105 


1.2522 x lO" 08 


2.1429 x 10+ 02 


1.5156 x 10+ 09 


9.4892 


13.9882 


0.0106 


7.1354 x 10~ 07 


2.1226 x 10+ 02 


2.6708 x 10+ 07 


9.5286 


13.8645 


0.0107 


2.5363 x 10~ 05 


2.1028 x 10+ 02 


7.5462 x 10+ 05 


9.5697 


13.7402 


0.0108 


5.3391 x lO" 04 


2.0833 x 10+ 02 


3.6009 x 10+ 04 


9.6129 


13.6151 


0.0109 


6.1795 x 10" 03 


2.0642 x 10+ 02 


3.1259 x 10+ 03 


9.6583 


13.4888 


0.0110 


3.5077 x lO" 02 


2.0455 x 10+ 02 


5.5342 x 10+ 02 


9.7062 


13.3610 


0.0111 


8.4175 x lO" 02 


2.0270 x 10+ 02 


2.3183 x 10+ 02 


9.7570 


13.2314 


0.0112 


9.5984 x lO" 02 


2.0089 x 10+ 02 


2.0443 x 10+ 02 


9.8112 


13.0996 


0.0113 


1.2079 x lO" 01 


1.9912 x 10+ 02 


1.6341 x 10+ 02 


9.8692 


12.9648 


0.0114 


1.3117 x 10~ 01 


1.9737 x 10+ 02 


1.5143 x 10+ 02 


9.9318 


12.8264 


0.0115 


1.4639 x lO" 01 


1.9565 x 10+ 02 


1.3662 x 10+ 02 


10.0000 


12.6834 


0.0116 


1.6190 x 10~ 01 


1.9397 x 10+ 02 


1.2446 x 10+ 02 


10.0752 


12.5344 


0.0117 


1.7538 x lO" 01 


1.9231 x 10+ 02 


1.1586 x 10+ 02 


10.1594 


12.3773 


0.0118 


1.8752 x lO" 01 


1.9068 x 10+ 02 


1.0938 x 10+ 02 


10.2559 


12.2088 


0.0119 


1.9940 x 10" 01 


1.8908 x 10+ 02 


1.0402 x 10+ 02 


10.3703 


12.0232 


0.0120 


2.1463 x 10~ 01 


1.8750 x 10+ 02 


9.7983 x 10+ 01 


10.5150 


11.8082 


0.0121 


2.2964 x lO" 01 


1.8595 x 10+ 02 


9.3439 x 10+ 01 


10.7287 


11.5252 



TABLE III: The computed values of TP in t, PE max , r, au v and a rtp for 22 different values of A when E — 185. 



